Thermally Induced Density Perturbations in the Inflation Era 
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The possibility of thermally induced initial density perturbations in inflationary cosmology is 
examined. The fluctuation dynamics of a scalar fleld plus thermal bath system during slow roll 
is described by a Langevin-like equation. Fluctuation-dissipation arguments show that for a wide 
parameter range within the standard inflation model, the thermal fluctuations of the scalar field can 
dominate its quantum fluctuations. The initial amplitude of density perturbations is found to lie in 
a range which is consistent with the recent observations of cosmic temperature fluctuations. 

PACS numbers: 98.80.Cq, 05.40. +j 

■ The standard scenario for structure formation in the 
' universe is based on inflationary cosmology. According 
to this model, quantum fluctuations of the scalar field 
during the expansion era were the perturbing seeds in 
^ . an initial, globally smooth universe. From this, large- 
scale structures then arose [||. This model predicts that 
the initial density perturbations should be Gaussian, and 
have a power-law spectrum with index n ~ 1. Based on 
I , this model for the initial perturbations plus assumptions 
^ ■ about dark matter, the formation of galaxies, clusters of 
■rj" galaxies and other objects have been extensively stud- 
ied 1^. Subsequent studies have shown that this model 
is consistent with observation, thus not theoretically im- 
plausible. Recent detection of temperature anisotropies 
in cosmic background radiation (GBR) by the Differen- 
, tial Microwave Radiometer (DMR) on the Gosmic Back- 
' ground Explorer's {COBE) satellite has given the first 
opportunity to directly probe the initial density pertur- 
bation. These results fit the scaling spectrum given by 
O ' the infiation model j|] . Thus it is one source of support 
^ , for the infiation model in describing the initial density 
^ ' perturbations. 

. . ' However, the COBE-DMR results do raise questions 

^ ' about the amplitude of the initial perturbations. Such 
k>( ' questions are not new, only further perpetuated by the 
COBE results. Before COBE it was already possible 
to determine the amplitude by fitting (or normalizing) 
the evolved perturbations with the observed clustering of 
galaxies on scales of say 8 h~-'^Mpc, where h is the Hub- 
ble constant in unit of 100 km s~^Mpc~^. This method 
suffers from three uncertainties arising from: 1) the evo- 
lution of the perturbations; 2) the assumptions of dark 
matter, and 3) the bias factor. On the other hand, the 
temperature fiuctuations of GBR on scales of superhori- 
zons size at the decoupling time, i.e. larger than about 
2°, directly provide information about the initial den- 
sity perturbation, independent of the above three factors. 
The COBE result of the GBR quadrupole amplitude has 
already been used to normalize the amplitude of the cos- 
mic density perturbations. Therefore, it is now a ques- 
tion of observational importance to explain the so found 
amplitude of the initial perturbations. 

The origin of the amplitude's magnitude is also rele- 
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vant because the "standard" scenario cannot explain it 
naturally. It has been known since developing the in- 
flation model, that the amplitude of the initial density 
perturbations given by quantum fluctuations of the infla- 
tionary scalar field is, at least for the "standard" model, 
of order 1, which is about 10* times larger than the re- 
quired value [|j . This is sometimes called the fluctuation 
problem. One can relax this conflict by sophisticated de- 
signs for the inflationary potential. However, such par- 
ticle physics designs counter the naturality philosophy 
of inflation. Moreover, in these models the amplitude is 
completely determined by unknown parameter (s) of the 
so designed potential. This gives almost no constraint 
to the possible or reasonable range for the amplitude. 
As such, it cannot predict what the initial value of the 
perturbation amplitude should be. 

This situation motivates a search other possible mech- 
anisms, which do not depend on such sophisticated de- 
sign. In this paper we show how thermal fluctuations 
during inflation may actually play the dominant role in 
producing the initial perturbations. It is conventionally 
believed that the components of particle-like matter, ei- 
ther relativistic or non-relativistic, are totally negligible 
during inflation. This is certainly true if we only consider 
the energy density, because inflation is by definition the 
epoch when the vacuum energy of the scalar field was 
the dominant component in the universe. However, this 
docs not imply that the initial perturbations must mainly 
arise from quantum fluctuation of the scalar held. All 
particle-like matter which existed before inflation would 
have been dispersed by inflation. Yet, particle-like matter 
will not completely vanish if one considers the processes 
of the scalar fleld dissipating into a thermal bath via its 
interaction with other flelds. 

The existence of a thermal component during infla- 
tion may not be exceptional and perhaps even inevitable. 
Roughly one can see this as follows. In order to maintain 
the (f) field close to its minimum at the onset of the infia- 
tion phase transition, thermal forces will generically be 
an important contributing source. Therefore, at least in 
the starting period of the phase transition, there is ther- 
mal contact between and all other fields with which it 
interacts. During the slow roll period of inflation, the ki- 
netic and potential energy of the (j) field is fairly constant, 
so that the interaction between the field with the other 
fields remains about the same as at the beginning. As 
such, there is no compelling reason to believe that the 
thermal component vanishes during inflation. 

We will show that, for a wide range of the parameters 
characterizing inflation, there can be a thermal compo- 
nent in the universe. With account for such a compo- 
nent, we flnd that the density perturbations can origi- 
nate mainly from thermal fluctuations. Accepting this 
mechanism, the fluctuation problem would be automat- 
ically avoided. In addition the allowable range for the 
amplitude would be consistent with the observed results. 

Let us consider the standard model of inflation, given 
by a scalar field (j) with Lagrangian density 
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L = -d^^df,^-V{(^)+L,nt , (1) 



where is the effective potential, and Lint describes 

the interaction of (p with all other fields. The classical 
equation of motion for in a de-Sitter universe is 

4> + 3H^ + T^<j)-e~^"'S7^(j) + V'{^) ^0 . (2) 

The friction term T^cji phenomenologically describes the 
decay of the </> field via the interaction Lagrangian Lint- 
In principle, the friction term T^(j) may not be reason- 
able for describing the energy transfer from 0-field parti- 
cle production during far out of equilibrium conditions. 
However, it is a proper approximation for describing the 
energy dissipated by the 4> field into a thermalized radi- 
ation bath. It should also be noted that could be a 
function of </>. Due to the lack of a detailed model for the 
decay of 0, we assume that is a constant independent 
of (j>. As we will show, the thermal fluctuations of <j) do 
not depend significantly on the details of F^. 

In the standard treatment of the inflation model, one 
assumes that the inflation era is divided into two regimes: 
1) slow roll and 2) reheating. In the former all inter- 
actions between the inflationary scalar field and other 
fields are typically neglected. These interactions have 
been considered only in the latter, in order to supply the 
mechanism to reheat the universe. This is equivalent to 
assuming that the friction term is only important in the 
reheating regime, but negligible during the slow-roll (in- 
flation) regime, i.e. 3H 3> F^. In this case, the slow- roll 
evolution in a spatially homogeneous universe is given by 
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(3) 



where = (87rG/3)p0 - (87r/3)AfV"^p;, P'P = '/'V2 + 
y(^) is energy of the fleld, rripi = is the 

Planck mass, and V{(j)) ~ V"(0) = M'^. Eq.(3) is 
valid when the potential V{(p) satisfles the following well 
known conditions for the inflationary potential , 

\V"{4>)\ « 2W(0)/m2;, 

V'\cP)m% « 487rVH(t>) ■ (4) 

Our flrst observation is that the condition F^ ^ 3H 
is not necessary for a slow roll solution. The coupling of 
the inflationary scalar fleld with other flelds can co-exist 
with the roll down solution. Consider the case when F^ 
is comparable to H. This implies that its decay prod- 
ucts will equilibrate quickly to some temperature T^. For 
explicitness in the treatment below, let us make the rea- 
sonable assumption that the decay products of the 0-fleld 
are relativistic matter. The additional equation needed 
to describe this relativistic component from the first law 
of thermodynamics is 

Pr + AHpr = F0(/)^ , (5) 
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where pr is the energy density of the thermal component. 
With account for this component, the slow roll equation 
(3) should be replaced by the set 
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Pr ^ 0. (7) 

and 

if' = — (P0 + Pr) . (8) 

Strictly speaking, when there is a thermal component in 
the universe, one should use a finite temperature effec- 
tive potential to replace the zero-temperature potential 
V{4)). As we will show below, this replacement will not 
be important for our purpose. 

Solution (7) implies that independent of the initial con- 
ditions for the thermal component, it will reach a steady 
state regime during inflation, i. c. the depiction of the 
radiation due to expansion will be balanced by its pro- 
duction due to friction. Prom eqs.(5) - (7), the constant 
energy density of the thermal component is found to be 

Pr - . (9) 

In the inflation epoch the kinetic energy of the (f) field, 
0^/2, is much less than its vacuum energy ~ ^{4')- 
Thus we have pr -C if 

< aAH , (10) 

where a > 1 is a model dependent arbitrariness. As such, 
in terms of energy density, the considered system during 
inflation is still dominated by the vacuum energy of the 
(f) field, with in particular the thermal component being 
negligible. Thus other aspects of the the inflationary sce- 
nario will remain the same as in the standard model. 

However, in terms of the system's temperature, the 
presence of a thermal component is not necessarily neg- 
ligible. The temperature, T^, of the thermal component 
is given by 

Tr ~ pI'^ ~ {mj,iWT^Y'^M^I^ , (11) 

where W = (fP' /2V{4)) is the ratio of the kinetic and po- 
tential energy of the (j) field- From eqs. (9), (10) and 
(11), it is easy to show that the temperature of the ther- 
mal component can be greater than the Hawking tem- 
perature, i.e. 

Tr>H (12) 
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All conditions (4), (10) and (13) can be simultaneously 
satisfied if 



V^'\^)ra-J « F'(0) « m-^V{^) (14) 

Since (M/nipi)'^ <C 1, this condition can be fulfilled. In 
fact, there is big room in the parameter space of the 
potential, in which the inflationary expansion of the uni- 
verse will still be dominated by the scalar field, but the 
temperature of the system will be determined by the ther- 
mal matter pr. On the other hand, because W 1 dur- 
ing the slow-roll regime, eq.(ll) implies that <^ M. 
Using A(/)^ as a generic inflationary potential, the lead- 
ing temperature effect is known to be XT^ [Q. Recall- 
ing that to have inflation requires A <C {M /mpi)^ means 
AT^ < M^/mli ~ H'^. Therefore, the influence of the 
finite temperature effective potential is insignificant as 
stated earlier. 

The quantum mechanical fluctuations of the (j) held 
during inflation are determined by the Hawking temper- 
ature H. Therefore, one can expect that eq.(12) is the 
condition under which thermal fluctuations will compete 
with quantum fluctuation. We will justify this point in 
the following. 

To calculate the fluctuations, the (f> field should be 
treated as a stochastic field. Therefore, eq.(6) should be 
interpreted as the ensemble averaged equation of motion. 
The essence of eq.(6) can easily be seen if it is rewritten 
as 

# ^ 1 dF[4>] 

dt 3H + T^ d(f> ' ^ ^ 

where F[(l)] = V{(j>). Eq.(15) is, in fact, a equation for 
the rate of change of the order parameter of a homo- 
geneous system with free energy F[(j)]. It describes the 
approach to equilibrium for the system. More generally 
in a study of fluctuations, one should not use the approx- 
imation of a spatially homogeneous universe, so that the 
spatial gradient term cxp(— 2i7t) in eq.(2) should not 
be ignored. The equation (15) for the rate of change is 
then modified to 

#(x,0 _ 1 ^f[0(x,t)] 

dt 3H + T^ 5(j) ' 

where the free energy is given by 

1 



-(e--W+n0) 



(17) 



Our present purpose is limited to comparing the ampli- 
tudes of thermal and quantum fluctuations. For this, it 
is enough only to calculate the 0-field fluctuations for 
the mode of wavelength equal to the horizon H^^ during 
inflation. As such the contribution from the spatial gra- 
dient term for only this mode is important in calculating 
the correlation function of the scalar field. 

It is known that rate equations like (17) cannot cor- 
rectly describe the approach to equilibrium during a 
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phase transition without also a noise term For 
instance, eq.(17) will only cause the order parameter 
to evolve towards local minima but not necessarily the 
global minimum. To ensure that the system approaches 
the global minimum, we must remember that actually the 
order parameter dynamics is not purely relaxational, but 
may exhibit fluctuations, arising from the microscopic 
degrees of freedom. These fluctuations can be modeled 
by introducing a noise term i] into eq.(17) as 

d(t) 1 5F , , 

-+rjit). (18) 



dt 3H + S( 

This is the Langevin equation for a system with one de- 
gree of freedom, in similar analogy to say a Brownian 
particle. A similar type of equation has been examined 
in 1^. However, their purpose was to statistically treat 
the quantum fluctuations of the scalar field, whereas ours 
is to treat external thermal forces. 

We had shown in ^ that during the eras of dissipa- 
tions, the dynamics of structure formation in the uni- 
verse should be described by a KPZ-equation which 
describes systems governed by non-linear effects plus 
stochastic fluctuations. Eq.(18) is a realization of this 
hypothesis, within the context of scalar field dynamics 
in the standard inflationary model. Although it is in all 
generality non-linear, in this paper we will not study the 
complete spectrum of the density perturbations, but only 
the amplitude of fluctuations with wavelength ~ H~^. 
For this we will not concentrate on the non-linear effects 
in eq.(18). 

The stochastic force rj in eq.(18) can be found from the 
fluctuation-dissipation theorem Q. If the temperature 
of the thermal bath is T^, the expectation values of 77 is 
given by 

m), = (19) 

and 

{v{tHt'))v = DS{t - t') . (20) 

The notation (...),, denotes averaging of 77 with respect 
to a Gaussian distribution. The variance D is given by 

1 Tr 

D^2 (21) 

where U — {4:Tt/3)H^^ is the volume with Hubble radius 

The fluctuations S(l) of the 4> field can be found from 
linearizing eq.(18). If we only consider the fluctuations 
d(j) crossing outside the horizon, i.e. with wavelength ^ 
H^^. The equation of Sip is 

^^ = -^^±ilM^^-f,. (22) 

According to condition (4), the term V {(/>) on the right 
hand side of eq.(22) can be neglected. From eq. (22), 
one finds for the correlation function of the fluctuations 
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(23) 



so that 



(('5'/>)')'7 - ^HTr . (24) 

This is our central result. Notice that it is independent 
of r^. This is expected since it is simply the variance 
of the (p field when coupled to a thermal bath, as im- 
plicit to the fluctuation-dissipation theorem. Neverthe- 
less, the importance of a sufficiently large decay term, 
as emphasized earlier, is to ensure appropriate dynam- 
ical conditions for rapid thermalization of the radiation 
bath on the scale of the expansion rate H. From eq.(24) 
one can conclude that the thermal fluctuations of the 
scalar field will be greater than its quantum fluctuations, 
{{5(t>)'^)QM ~ -ff^/27r, when condition (12) holds. 

Since the kinematics and dynamics of inflation arc the 
same here as in the standard model, the initial pertur- 
bations will still have a power-law spectrum with index 
n \. Also the amplitude 5p/ p, when it crosses back in- 
side the horizon, can be calculated by the gauge invariant 
amplitude Q = 5p(j,/{p + p) during the time of inflation. 
For quantum fluctuations it is known that C is of order 
1. This follows simply because for a field in its ground 
state, the mean quantum fluctuation of its energy density 
is of the same order as its mean kinetic energy density. 
However, for thermal fluctuations, the mean kinetic en- 
ergy density will be greater than its fluctuation, so one 
can except that the quantity C, should be much less than 
1. 

The energy fluctuations caused by 5<j) is 5p = 5(j)V'{4i), 

and p + V = P<t> + V<t> + Pr + Pr = 4>^ + (4/3)pr-- Therefore 
from eqs.(6), (9), (11) and (24), we have 

V ^ ^l(E±\^^ (eX'^ (25) 

P <P^ + {4/3)pr A.\t,h) \TrJ 

Because Tr < M, eq.(25) shows that the amplitude of 
the initial perturbations, Sp/p, should be in the following 
range. 



rripij 



< - « 1 • (26) 



Therefore, the amplitude of the initial perturbations is 
mainly limited by the ratio M/nipi, i.e. the energy scale 
of inflation. If we take M ~ lO^^Geu, the possible range 
for the amplitude 6p/p should be in about the middle of 
the range (10^*^ - 1){T^/Hy/'^. This result is consistent 
with the observed amplitude Sp/p ^ 10^'*. 

As an additional outcome of this treatment, eq. (26) 
places an upper limit to the energy scale of inflation of 
M < m,pi{Sp/ p)'^^'^ ^ m,pi ■ 10~^, above which thermally 
induced fluctuations would be inconsistent with the ob- 
served density perturbations. Therefore, one can also 
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conclude that for thermally caused initial perturbations, 
inflation should not occur earlier than about TOp//10^. 

In addition to the amplitude fluctuation for the scalar 
mode, which was treated in this paper, there is also an 
amplitude fluctuation for the tensor mode . Since this 
involves weakly interacting gravitons, a thermal mecha- 
nism for inducing these fluctuations seems less likely to 
the standard treatment which considers quantum fluctu- 
ations. 
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